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^ ' Abstract 
> 

t^. | We investigate the dynamical 3-form flux compactifications and their implications for 

| cosmology in our brane world in the six-dimensional Nishino-Salam-Sezgin model, which is 

, usually referred to as the Salam-Sezgin model. We take the background of the 3-form field 

acting on the internal space and timelike dimensions without the U(l) gauge field strength. 
The first class of solutions we discuss is the dynamical generalization of the static solutions 
obtained recently. In this class, we find that the time evolution is restricted only in the 
^ | azimuthal dimension of the internal space and not in the ordinary three-dimensional one, 

which does not give a cosmological evolution. The second class of solutions is obtained by 
exchanging the roles of the radial coordinate and the time coordinate from the assumptions 
in the first class. At the center of the internal space, there is a conical singularity which 
may be interpreted as our 3-brane world, but it is difficult to realize a warped structure 
in the direction of the ordinary 3-space. Except for the oscillating solutions, a dynamical 
evolution leads to expanding or contracting ordinary 3-space, depending on the choice 
of time-direction. Furthermore, in the expanding solutions, there are decelerating and 
accelerating ones: In the latter solution, the evolution is sustained for a finite proper time 
and the scale factor of the 3-space diverges. 
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1 Introduction 



Brane world models have attracted particular interests in recent years. In the cosmological 
aspects, six- dimensional brane world models may be useful for a resolution of the cosmological 
constant problem [1] (see, however, [2]). A six-dimensional model is also recognized as an 
important playground to study cosmology and gravity with stabilized extra dimensions by fluxes 
of antisymmetric tensor fields. In the simplest realization of the flux compactifications in six 
dimensions, the internal space has the shape of a rugby ball [1,2], where codimension-two branes 
are placed at the poles. The warped generalizations of the rugby ball solutions were reported in 
the Nishino-Salam-Sezgin (NSS) six-dimensional supergravity [3] and in the Einstein-Maxwell 
theory [4]. See [5] for the NSS model. Note that this model is usually referred to as the Salam- 
Sezgin model especially when compactification is discussed. However, the compactification on 
S 2 , which was found by Salam and Sezgin, was constructed in the theory first proposed by 
Nishino and Sezgin. Therefore, it is more suitable for us to refer to the model as the NSS 
model. In Ref. [6], a study of the cosmology was presented on the analogy of the classical 
mechanics. It is also recognized that a 3-brane in six or higher dimensions faces a problem 
about the localization of the matter, because such a brane generically produces a curvature 
singularity. A way to circumvent this problem is to introduce a thickness of the brane, as 
proposed in [7, 8]. 

In order to explore the cosmology, one useful approach is to use the dynamical solutions. In 
this direction, exact dynamical solutions were obtained in [9, 10, 11, 12]. The field components 
of the bosonic part of the NSS model are the gravity, the ?7(l)-gauge field, Kalb-Ramond 2-form 
(i.e., 3-form field strength) and the dilaton. In Refs. [9, 10], warped time dependent solutions 
with the [/(l)-field strength, without the Kalb-Ramond 3-form field strength in the NSS model, 
were explored. The solutions have all six dimensions evolving in time. There are three types of 
the exact time dependent solutions: One of them is the scaling solution, in which both the scale 
factors of the external and internal spaces have the linear dependence on the proper time. The 
other two solutions have more complicated time dependence. In the earlier times, these solutions 
can be seen as the generalizations of the solutions obtained in Ref. [13], but in the later times 
both approach the above scaling solution. The scaling solution is the generalization of the static 
solution obtained in [3] and becomes an attractor. In this sense, the static brane solution in 
the 2-form background is always unstable. An extension of the time dependent 2-form solution 
in the NSS model to more general Einstein-Maxwell-dilaton models was considered in Ref. [11]. 
The exact wave solution was also obtained in [12]. 

The purpose of this work is to present the other classes of time dependent solutions in the 
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NSS model in the presence of the nonvanishing 3-form field strength [but in the absence of 
the U(l) field strength], and discuss their implications for the brane world cosmology. The 
first class of solutions is the dynamical generalization of the static solutions obtained recently 
in Ref. [14], where the internal space has a torus topology with two cusps, along the line of 
the similar generalization of brane solutions [15]. We will find that our solutions are stable, 
in contrast to the 2-form solutions. The second class is obtained by exchanging the roles of 
the radial coordinate and the time coordinate in the first class. We also discuss the dynamical 
generalization of the black 1-brane solution in the NSS model. 

The paper is constructed as follows. In Sec. 2, we briefly review the NSS model. In Sec. 3, 
the time dependent generalization of the solutions and their application in constructing the 
brane world model are discussed. In Sec. 4, the S-brane-like solution, obtained by exchanging 
the roles of space and time, and their cosmological properties are discussed. The final section is 
devoted to give a brief summary and conclusion. In the Appendix, we discuss the black 1-brane 
solution. 



2 Nishino-Salam-Sezgin model 



The action of the bosonic part of the NSS supergravity [5] is given by 
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where <fi is the dilaton, H = dB is the field strength for the Kalb-Ramond field B and F = dA is 
the field strength for the U(l) gauge field A. The parameters g and k are gauge and gravitational 
coupling constants, respectively. Variation of the action (2.1) gives the field equations: 
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Most of the preceding works have discussed the U(l) field acting on the internal space di- 
mensions, F 7^ 0, with the vanishing 3-form H = 0. Instead, in this work, we consider the 
nonvanishing 3-form field acting on the internal space and timelike dimensions H ^ 0, with the 
vanishing U(l) field F = 0. 



2 



3 Generalized warped compactifications by the 3-form 

3.1 Ansatz and equations 

In this section, we take the following metric ansatz [15]: 

3 

ds l = - e 2u ^dt 2 + e 2u ^ ^(rfx 4 ) 2 + e 2v ^dy 2 + e 2w ^d6 2 : (3.1) 

i=i 

where the coordinates t and x l describe our four-dimensional worldvolume, and the remaining 
y and 9 are transverse to it. The metric components u , ui, v, w and the dilaton <p are assumed 
to be functions of both t and y. For convenience of our derivation, we also define 

U = u + 3ui — v + w. (3.2) 

We assume that for the 3-form field strength there is only the y dependence and also the U(l) 
field strength vanishes: 

H = E'(t, y) dt A dy A d9, F = 0. (3.3) 

Throughout the main text in this paper, the dot and prime denote derivatives with respect to 
t and y. The Einstein equations (2.2) become 



(«o-«)(ug + U'u'q) -U + u -2v + u (U -ii + 2v) - 3ii\ - v 2 - w 2 



Su'i — w' + u' (3ui + w) — ?>u\{ui — v) + w'(v — w) — cfxp', (3.5) 
e2(«i-«)( u // + u' u '^ _ e - 2 («o-«i)[ iil + ^(U - 2ii + 2v)} 

2 2 
= _^_ e -2<t>-2v-2w £/2 _ 9_ e <t>+2uu ^3 g\ 

U" + v" + u' 2 + 3u[ 2 + w' 2 - v'(U' + v') - e- 2{u °- v) [v + v(U - 2ii + 2v)] 

= _0' 2 + ^g-20-2,0-2^/2 _ (3 ?) 

2 K 

e 2(— + u' w ') _ e 2 (—«o)^ + ^(f/ _ 2u + 2v)] 

2 2 
_ ^ e -2^2« -2i>^2 _ 9_ e 4>+2w^ ^3 



The dilaton Eq. (2.3) and the equation for the 3-form field (2.4) give 
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( e tyV _ f e U-2uo+2v^y _ K 2 e -n +3 Ul -^-2^/2 _ ^9^ = ^ ,3 
( e U-2u -2 W -2<t, E y = Qj (3 
( e U-2u -2 W -2t E y = Q (3.11) 



The latter Eqs. (3.10) and (3.11) give 

e U-2uo-2w-2<j> E , = ^ ( 3 _ 12 ) 

where c is an integration constant. Substituting (3.12) into Eqs. (3.4), (3.6), (3.8) and (3.9), we 
obtain 

e u u' - y cE) ' = e u+2{v - Uo) [U-u + 2v- ii (U -ii + 2v) + 3u\ + v 2 + w 2 + 2 ] 

-^ +2v+u , (3.13) 

K 

2 2 
e U v! 1 + —cE\ = (u ie U+2(v-u )y _ 9n2v + U ^ ,3^ 
2 / K 

VV-^cflY = (uie" 7 ^)" - 4e^ +c/ , (3.15) 
2 / /t 

?<!>'- k*cE)' = (0 e ^+2(-«o)y + V e ^+cr_ (316) 

/ K 

Let us assume that the solutions satisfy the following equations: 
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V e< H2„ +f / _ 2£ , (3.20) 

ft 

where li ii = 0, 1, w, <fi) are constants corresponding to the time derivative parts in Eqs. (3.13)- 
(3.16). Relaxing these separability conditions is a very interesting issue, but we find it difficult 
to do so. As in Ref. [15], we also assume that U is independent of y. In fact, it is known that 
some nontrivial time dependent solutions can be obtained even with this restriction. Then using 
v' = u'q + 3u[ + w' obtained from U' = and Eqs. (3.17)— (3.19), we find 

C„2 



(2 V + <j>)" = -JL e 2v+<t> + 2(£ + Mi U)e~ u . (3.21) 

It is hard to find a solution unless we assume 

£ + 3h + £ w - 1+ = 0. (3.22) 
When this is obeyed, the solution to Eq. (3.21) is 

e _(2^ )= V co8h 2 [/i(y-yi)], (3-23) 
K Ji 
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where fi and y\ are constants. Although relaxing the condition Eq. (3.22) would be also an 
interesting issue, we leave it for future study. 
It follows from Eqs. (3.17)-(3.20) that 

2 f 

u' = ^-ce- u E-^t a vh[f 1 {y-y 1 )]-G +£ e- u y, (3.24) 

2 f 

u[ = -^ce- u E-^tanh[f 1 (y-y 1 )]-G 1 + £ 1 e- u y, (3.25) 

w' = ^ce- u E-^taDh\f 1 (y-yi)]-G w + e vj e- v y, (3.26) 

0' = K 2 ce -u E + L tanh lf 1 (y- yi )] + 2G (j) -2^e- u y, (3.27) 
where we have chosen Gi to be constant for simplicity. We then find from (3.23) and (3.27) 

v 1 = -^ce~ u E - ^ tanh[A(y - y x )\ - G + ^e~ u y. (3.28) 
The condition U' — imposes 

G + 3d + G w - G^ = 0. (3.29) 

For convenience, we define 

Y = 2K 2 ce~ u E. (3.30) 
By substituting the above relations into Eq. (3.7), we find 



l(Y 2 -Y')-f 1+ e- u (l 



° 2 

+ 3A 2 + Al + 3A 2 + (3A 1 + A w -A <t> ) 2 + Y(3A 1 + A <t> ), (3.31) 

where we have assumed (ye u ~ 2uo+2v ) = £ v (ji v j s a constant), and have defined 

Ai = d - £ie- u y, (i = l,w,(f>). (3.32) 

At this stage, except for the cases where ii = 0, it is still not easy to find analytically tractable 
solutions. Thus, in this paper, we focus on the case ii = 0, leaving study of more general cases 
for the future. Then, without loss of generality, we may set 3Gi + G^ = 0. Now Eq. (3.31) 
reduces to 

Y 2 — Y' — f 2 = 0, (3.33) 

where we have defined 

39 

f 2 = 2f 2 - 4(G, - G w ) 2 - 6 -Gl (3.34) 
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f 2 can take either positive, zero and negative value and the solutions for Y(y) can be classified 
into the three types, depending on the sign of /J. We will give the solutions shortly. Once Y(y) 
is obtained, by integrating Eqs. (3.24)-(3.28), we can find the metric components and dilaton 

as 

n = -lp(y)-l\ncosh[f 1 (y-y 1 )]-(2G 4> -G w )y + h (t), (3.35) 
u x = ^P(y)-^\ncosh[f 1 (y-y 1 )} + ^G (t> y + h 1 (t), (3.36) 
v = ^P(y)-^lncosh[f 1 (y-y 1 )]-G <f> y + h v (t), (3.37) 
w = -^P{y)-^]ncosh[f 1 {y-y 1 )]-G w y + h w {t), (3.38) 

= -ip( y ) + ilncosh[/ 1 ( U - Wl )] + 2G^ + ^(t), (3.39) 

where the function P(y) depends on the sign of /f ; hence, the type of solutions, and hi (i = 
0, 1, v, w, <p) could be functions of t. When all hi vanish, the solution becomes static. 
Now the three types of solutions are given as follows: 

The sinh solutions: For a positive /|, the solution to Eq. (3.33) is given by 

Y(y) = -f 2 coth[f 2 (y-y 2 )]. (3.40) 

After integration, we obtain 

P(y) = \n\smh[f 2 (y-y 2 )}\. (3.41) 

The sin solutions: For a negative /|, the solution to Eq. (3.33) is given by 

Y(y) = -\f 2 \cot[\f 2 \(y-y 2 )}. (3.42) 

After integration, we obtain 

P(y) = ln\sm[\f 2 \(y-y 2 )]\. (3.43) 

The linear solutions: For /| = 0, the solution to Eq. (3.33) is given by 

Y(y) = — . (3.44) 

y-yi 

After integration, we obtain 

P(y) = \n\y-y 2 \. (3.45) 

These give the dynamical generalization of the static solutions found in [14]. The correspon- 
dence is given by the following reparameterizations of f\ -H- A 2 , f 2 <H> Ai and q -h- c. Note, 
however, that Eq. (3.34) does not completely agree with Eq. (2.10) of [14], in some factors. 
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3.2 General time dependent solutions 

Now our task is to solve the remaining time derivative parts of (3.14)-(3.16) and (3.7). The 
remaining time dependent parts of these equations give 

hi = (h - 3hi - h v - h w )hi, (3.46) 

where i = l,v,w,(f>. They lead to 

^ e _ ho+3hl+hv+hw = ^ ^ 34? ^ 

where all c$ are constants. The off-diagonal components of the Einstein equation (3.5) imposes 
the constraints among the time dependent functions 

h v = -3hi, h < p = Gh 1 , (llG<f> - 3G w )h 1 = (G w - G^)h w . (3.48) 

It is straightforward to confirm that except for the case of G w = G^ all hi cannot have 
any nontrivial time dependence. To show this, for the moment let us set h (t) = 0, which is 
always possible by an appropriate rescaling for the time coordinate. From Eq. (3.48), we have 
c v = — 3ci and = 6ci. 

In the case of 11G^ ^ 3G W , Eq. (3.47) tells us that 

■ G w G(fr • Gm Gfh 

fl\ = . — — - — r— , tl v — — O- 



^11G> — 3G w )t (11(j0 — 3G w )t 



It follows from Eq. (3.13) that 



U " (llG.-SGU 2 * 2 ' ^' 5Uj 



Thus, the consistency requires G w = G^. In the case of HG^ = 3G W , h w = 0. We may set 
h w = and then 3h\-\-h v = const. From Eq. (3.47), we have hi = Ct, h v = —3Ct and h^ = QCt, 
where C is a constant. We then find from Eq. (3.13) 

C 2 = 0, (3.51) 

which leads to a trivial solution. Thus, G w = G$ must be satisfied to have nontrivial solutions. 
For G w = (j^, 

fl = VI - fGl (3.52) 
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Then, h v = — h\ = but /i^, ^ 0. Thus we obtain 

hi = ki, h v = k v , h<f, = k^, (3.53) 

where hi, k v and k^ are integration constants. Therefore, the ordinary 3-space metric and the 
dilaton do not depend on time. From Eq. (3.23), we obtain 

^ = -2^ + m(-^)Y (3.54) 

h and h w are still undetermined except for the relation (3.47), which gives (e hw )' = c w e kv+3kl e h ° . 
The two-dimensional metric in the t and 9 directions now becomes 

ds * = e 2F(y)(_ e 2h 0dt 2 + e 2h wde 2^ 
p 2F(y) / 2 



= jL,^ i- dT2 + ey (t -+ 3fa) rW) , (3^55) 

where F(y) is the common y dependent part and T = e hw ^ is the proper time. Thus, the 
time dependence only appears in the 9 direction and not in the ordinary 3 directions, which 
cannot present any cosmological evolution. These results imply that the solution is dynamically 
unstable for the evolution along the 9 direction, but stable for the others. In terms of the four- 
dimensional effective theory, the effective potential would contain one flat direction associated 
with the motion in the 9 direction. 

One might guess that more generalized time dependent solutions could be obtained by 
allowing for the time variations of and Gj. However, in such cases, the time dependence of 
the Einstein and dilaton equations cannot be separated out from the dependence on y. Thus, 
it seems to be impossible to find consistent time dependent solutions. 



3.3 Comparison with 2-form flux compactification 

It would be interesting to compare our result with the dynamical solutions with the 2-form 
background in the NSS model. In the case of the U(l) background, there is the so-called scaling 
solution [9, 10], which can be the attractor. In this solution, both the internal and external 
dimensions have the time dependence linear in time and evolve uniformly. It means that the 
static solutions in the U(l) background are not stable. Thus, the scaling symmetry always 
requires some additional mechanism to stabilize the extra dimensions, for example, as discussed 
in [16, 17]. 
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In our 3-form background, the analyses in the previous subsection indicate that there may 
be no time dependent generalizations of the static solutions in the case of G w ^ G$. For the 
case G w = G^, the time evolution is allowed only in the 9 direction and is forbidden into the 
other directions. Thus in either case, this seems to imply that the 3-form compactification is 
stable in comparison with the case of £7(l)-background. Here one remark is in order. In deriving 
the solutions, we made some assumptions on the time dependence and relaxation of them might 
lead to additional possible evolutions of space. At the moment, we do not have the complete 
proof on the stability of such general solutions. 

3.4 Application to brane world model 

3.4.1 Singularity 

Reference [14] investigated the construction of the brane world by using the static solution in the 
background of the 3-form field strength. A curvature singularity exists at y = y 2 , irrespective of 
the type of solutions. Note that this is a real curvature singularity, which is more severe than a 
conical one and may not be identified as our 3-brane world. To see this, we define y :— y — y 2 , 
and the approximate spacetime metric near the singularity y = is given by 

ds* ~ cy/W + - + c^/ 2 jy x * f , (3 . 56) 

where {% = y, 6, t, x) are unimportant coefficients. This metric leads to a divergent scalar 
curvature R ~ y~ 5 ^ 2 near y = 0. Thus, as we will discuss in the next subsection, to construct 
a brane world, the singularity is first removed from the original spacetime in [14]. Then, an 
identical copy of the remaining piece of the spacetime is attached to the opposite side across the 
codimension-one boundary at y = y 2 + e, where e > 0, which wraps the axis of the rotational 
symmetry of the internal space. This boundary may be identified as our brane world. There 
are jumps of the physical quantities across the boundary. 

3.4.2 Construction of the brane world 

We now construct the brane world by employing the dynamical solution in the 3-form back- 
ground. First, the original spacetime is cut at y = yo > y\{> y 2 ) and the singular part of 
£/2 < y < yo is removed. Then, an identical copy of the remaining piece is glued with the 
original one, at the codimension-one boundary y — y$. 
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The induced metric on the boundary is given by 

dsl = e 2w{t ' yo) d6 2 - e 2uo(t > yo) dt 2 + e 2ui ^ yo) ^(d 
The extrinsic curvature to the boundary is given by 



.r 



i\2 



(3.57) 



Kab ~ 2e" ( 9ab ^ 



(3.58) 



</0 



where e = +1 or — 1 denotes the direction of the normal vector toward increasing or decreasing 
y, respectively. In our case, we choose e = +1. Now the extrinsic curvature is discontinuous 
across the boundary and its jump is determined by the Israel junction condition: 



9ab 



0. 



K. 



ab 



K, 



ab 



9abK 



(3.59) 



where [A] denotes the jump of a physical quantity A across y = yo, and S a b is the stress- 
energy tensor of the matter localized on the boundary. This co dimension-one boundary can be 
identified as our brane world. 

For our background metric, we obtain 



K% = -e~ v (u' + 3u[) , K\ = -e~ v (w' + 3^) , X -K\ 



S a b also can be decomposed into the components S^ = — p, S e e 
junction condition at y = yo is given by 

e ~ v (u' + 3%) 

'e^iw' + 3u' 1 ))\ v=yo 



y=yo 



K 

~2 



-P, 



e~ v (w' + u + 2u\) 



They can be rewritten as 



e- v (u[-w') 
e-"(ui + 3ui) 



\y=yo 



K 



2 y 



\y=yo 



\y=yo 



-y(p+p): 



K 

\y=vo = 



r v (u' + w' + 2u[). 

(3.60) 

and {1/3) SI = p. The 

(3.61) 
(3.62) 
(3.63) 

(3.64) 
(3.65) 
(3.66) 



Our discussion in this subsection can be applied to the sink and linear solutions discussed 
in the previous section, and the form of P(y) in Eqs. (3.35)-(3.39) is not specified. It is now 
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clear that the zero-thickness limit of yo — > 1/2 is not well behaved, since the left-hand side of 
the junction Eqs. (3.64)-(3.66) becomes singular. Note that the sin solution contains another 
curvature singularity at y = 7r / 1 y*2 j + 2/2 an d is excluded from our consideration. In addition, 
for the sink or linear solutions, the finiteness of bulk volume is not ensured. Thus, following 
Ref. [14], here we put the second brane world at some y = L > yo. Note that the second brane 
takes e = — 1. The junction condition at y = L is given in a similar way with the opposite sign 
of the right-hand side of Eqs. (3.61)-(3.63). 

In order to obtain a dynamical solution, we must choose G w = G^. Also, we may choose the 
integration constants as 

fci = 0, k v = 0, ^ = -21n(^fY (3.67) 

and thus {e hw ^j = c~ 1 e h °. By introducing the proper coordinate T = c w e hw , after some compu- 
tations, we obtain 

u l -w= l -P{y) + ^G <t> y- In T, (3.68) 

ux -u = -P{y) + ^y, (3.69) 

3«i + u = ^P(y) - In cosh [h{y - Vl )]. (3.70) 
The junction conditions are reduced to 

{e- V {l P '(yo) + ^)} = y (p e -p), (3.71) 

{ e ~ V (l P '(yo) + \g*)} = ~(P + P), (3-72) 

2 

{e- v (p'{y ) - fi tanh [frfo - Vl )]) } = \pe- (3.73) 

Noting that v is not an explicit function of T, the left-hand side of the junction equations become 
static, and p, p and pg remain constant. In the static case without the condition G w = G$, 
generically the three quantities p, p and pg satisfy the different equations of state from the 
tension. Now, however, from Eq. (3.71) and (3.72), the inclusion of the time dependence, i.e., 
the condition G w = G^, induces the relation p = —pg. 

Furthermore, the requirement of the recovery of the Lorentz symmetry at the brane y = yo, 
i.e., u[ = u' Q , together with Eq. (3.69) leads to 

3P'(y ) = -8GV (3.74) 

By combining this with the previous relations, u' = u[ = w' is obtained. This relation suggests 
that p = —p = —pg. Thus, we are lead to the conclusion that the brane with the ordinary 
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four- dimensional spacetime with the Lorentz symmetry is supported only by the tension. Note 
that this conclusion is derived with some assumptions concerning time dependence, and without 
them we might be able to construct the brane world with the Lorentz symmetry, supported by 
some matter other than the tension. 

Before closing this section, following Ref. [14], the four- dimensional effective theory on the 
brane world should be discussed. In the standard warped solutions, the four- dimensional Planck 
mass is obtained by integrating over the two-internal space. However, in our case, in general, 
the warp function in the timelike direction is different from that of the ordinary 3-space. Then, 
the six- dimensional Einstein-Hilbert term now reduces to 



Ml / (TxV^gR 



2(2tt)M 6 4 ! d 4 x^/^g~w( I ' dye^ U0+ ' iu ' +v+w g {A)tt R^ + f dye Uo+Ul+v+w g m R$ 
J Jyo Jyo 



(3.75) 

where M 6 := (7^2) l ^ is the six-dimensional Planck mass. Note that the factor 2 in front of the 
integration come from the Z 2 -symmetry, and (2ir) is from the integration over the azimuthal 
direction. To define the unique four- dimensional effective Planck mass M p , we need to require 
that 

2(2vr)M 6 4 / dye- U0+3ui+v+w = 2(2vr)M 6 4 / dye Uo+Ul+v+w =: M 2 , (3.76) 

Jyn Jyn 



or more explicitly 



dy , rr = / dy * , -=-. (3.77) 



cosh 2 [fx(y - yi)] J yo cosh 2 [} x (y - y ± )] 



In the case of ^ G w there is no time dependent solution. The cutoff parameter y appears 
in the definition of M 2 . In other words, the effect of the brane thickness is renormalized into 
the four- dimensional Planck mass. In the case of = G w , e w oc T and hence M 2 oc T. In this 
case, we should move to the Einstein frame. Then, the conformal transformation to the Einstein 
frame, gffi = Tgff, leads to the cosmic expansion of r 1 / 3 , where r is the cosmic proper time 
defined in the Einstein frame. However, this power in the expansion law is not sufficient for 
obtaining realistic cosmology. Realistic cosmological evolutions may be obtained by considering 
the time dependent matter on the brane, through the induced motion of the brane in the bulk. 
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4 S-brane-like solutions 



4.1 Field equation 

Assuming the same metric ansatz (3.1) as in the previous section, we investigate another class of 
the time dependent solutions by exchanging the roles of the t-coordinate with the y coordinate 
as well as those of uq with v. The way of constructions of solutions are very similar to the 
case of S-branes, see e.g., [18, 19, 20, 21]. We also assume that the 3-form field strength is the 
function of time 

H = E(t, y) dt Ady A d6. (4.1) 
The equations of motion for the flux are given by 



e 



U-2u -2w 



AJ-2u -2w- 



- 2 <t>E)' = 0, (4.2) 
- 2 <t>E)- = 0. (4.3) 



Equations (4.2) and (4.3) give 



where c is a constant and 



e 



$>-2v-2w- 



2<t> E = c, (4.4) 



$ := U + 2v-2u . 

In this section, we assume <3> = is independent of t, $ = 0. Employing (4.4), the Einstein and 
dilaton equations are now given by 

(eV)' = e*(u -u 2 + 3ul + v 2 + w 2 + 2 ) + \k 2 cE - 9 l e ^ +u , (4.5) 

2 K 

(eV)' = (^ e * " \k 2 cE)- - g l e ^ +u , (4.6) 

2 K z 



e 



u (v" + U" - v' 2 - v'U' + 3m/ 2 + w' 2 + u ' 2 + 0' 2 ) 



= fa* + 1 -k 2 cE)- - ^e^ +c/ , (4.7) 
2 K 

(e u w>Y = (we* + -k 2 cE)- - ^ +2v+u , (4.8) 

v ' 2 K 

(e u <j/)' = (0e* + K 2 cE\ + ^ e ^ 2v+u . (4.9) 
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Let us look for the solutions with the following assumptions that 

= e * (n - ul + 3u 2 + v 2 + w 2 + 2 ) + \k 2 cE - 9 l^ +2v+u , (4.10) 



K 

,2 



= fee* - -k 2 cE)- - ^ +2v+u , (4.11) 

2 K 

= (i) e * + I K 2 C £)- - 4e^ +c/ , (4.12) 



« 2 



= (ti;e $ + ^ 2 cE)- - ^- e ^ +c/ , (4.13) 
2 fir 



= (0e* + K 2 cE)- + ^e^ +c/ . (4.14) 

K 

By combining Eq. (4.11)-(4.14), with $ = and hence Uq = v + ?>ii\ + w, we find 

(2 Mo + 0)" = ^e 2 ^. (4.15) 

For Eq. (4.15), there are three classes of solutions, namely the sink, sin and linear solutions. 
We discuss these solutions in order. 



The sinh solutions: A solution to Eq. (4.15) is given by 

e -( 2uo +0) = V sinh 2 [/i(t _ (4 . 16) 



where f\ and ti are constants. Assuming that the right-hand sides of Eqs. (4.5)-(4.9) vanish, 
we obtain 

2 f 

v = -^ce-*E- J -±coth{h{t-h)\-G v , (4.17) 

2 f 

iii = ^-ce^E -^coth[f 1 (t-t 1 )]-G 1 , (4.18) 

2 / 

w = -^-ce^E -^coth[f 1 (t-t 1 )}-G w , (4.19) 
= _ K 2 ce -<f£ + A co th[/i(t-ti)]+2G^, (4.20) 

where Gj (z = i>, 1, w, <fr) are also constants. Because of the assumption that $ = 0, we obtain 

u = Y ce ~" E ~ x coth ^ - - ( 4 - 21 ) 

and 

3d - G + G„ + G w = 0. (4.22) 
11 



Substituting these into Eq. (4.5), we find 

^(Y 2 + Y)-(3G 1 + G (j) )Y 

—fx + ( — 3Gi + — Gw) + 3G\ + G 2 W + = 0, 
where Y is defined as in the previous section: 

Y = 2K 2 cEe-*. 

Without loss of generality, we may impose the further condition that 



3d + Ga = 0. 



Then Eq. (4.23) reduces to 



where we defined 



Y 2 + Y + f. 



2=0, 



39 

/I = -2A 2 + 4(G^ - G w ) 2 + y G}. 

Note that /f can be either positive, negative or zero. 

For a positive, negative and vanishing /|, the solution to Eq. (4.26) is given by 

^ = h cot[/ 2 (t - t 2 )], |/ 2 | coth[|/ 2 |(t - t 2 )], 



t-t 2 

respectively. The corresponding solutions with the possible y dependence are given 
"i 



-Q(t)-|ln 



smhlMt - h)} +-G 4> t + h 1 (y), 



w 



-~Q(t) - i In | sinh[A(t - ij)] | - (2G^ - G w )t + h v (y), 
-~Q(t)--ln sinhf/i^-^)] -G w t + h w (y), 



sinh[/i(t-ti)] + 2G t + /i (y), 
where for a positive, negative or vanishing /|, Q(t) is defined by 



Q(t) :=ln sin[/ 2 (t-t 2 )] 



In 



sinh[|/ 2 |(t-t 2 )] 



In It - t 



2b 



respectively. 
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The sin solutions: Another solution to Eq. (4.15) is the oscillatory type: 



-(2u +</>) 



Then, defining the new parameter 

fl = 2fl + 4(G, - G w f + 



, 32 2 



(4.35) 



(4.36) 



which is always positive for a positive , we get the corresponding solutions with the possible 
y dependence: 



u 

Ul 

v 
w 



I In 

4 

iln 
4 

-i In 

4 

> 

-iln 

2 



sin[/ 2 (t-t 2 )] 
sin[/ 2 (t - t 2 )] 
sin[/ 2 (t - 1 2 )] 
sin[/ 2 (t - 1 2 )] 
sin[/ 2 (t - i 2 )] 



5 In 
4 

iln 

4 

1 



sinLMt-tx)] -G^ + h (y), 



sin[/i(i-ti)] + + 



In 



In 



iln 

2 



sin[/!(t - t x )] - (2G^ - Gjt + ^(j/), 
sin[/i(t - *i)] - G w t + h w (y), 
smiMt-h)] + 2G <j) t + h (j) (y). 



(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 



The linear solutions: The last solution to Eq. (4.15) takes the form 

4g* 



-(2« +</>) 



K 



(t-hy. 



Then, by defining 



/I — 4(G 9 i — G h 



32 



Gz 



(4.42) 



(4.43) 



which is always positive, we get the corresponding solutions with the possible y dependence: 



"o 

V 

w 



iln 

4 

iln 

4 

-iln 

4 

-iln 

2 



sin[/ 2 (t 


-**)] 


--In 


sin[f 2 (t 


-**)] 


-llnlt-txl 



where 



-G^t + /ii(?/), 



sin[/ 2 (t - t 2 )] - - In |t - ti| - (2G^ - Gjt + /i„(y), 
sin[/ 2 (t - t 2 )] - - In |t - ti| - G w t + h w (y), 
sin[/ 2 (t - t 2 )] | + - In |t - tx| + 2G t + 



f 2 = vi + m - g w ? + f g\. 



(4.44) 
(4.45) 
(4.46) 
(4.47) 
(4.48) 

(4.49) 



In all the three types of solutions, setting all hi = leads to purely time dependent solutions. 
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4.2 Generalization to the y dependent cases 

Following the similar arguments as done in the previous section, we find that the y dependence 
can be included only for the choice of G<j> = G w . In any case of the sinh, sin and linear solutions 
shown in the previous subsection, the y dependent functions satisfy 

h = k , h 4> = k h h 1 =k 1 , (e h ™)' = c w e- k °-^e h % (4.50) 

where ko, k v and k\ are integration constants satisfying 

k^ + 2^ = -2 In ( ) . (4.51) 



Thus, even taking into account the y dependence, the spacetime structure in the ordinary four 
dimensions does not contain a warped structure. The internal space metric can be written as 

ds\ = F(t) 2 (e 2h *dy 2 + e 2h "d6 2 ) 



C w 



p 2(feo+3fei) p( f \2 / 2 

' dR 2 + -^,R 2 de 2 ) , (4.52) 



P 2 \ p2(fco+3fci) 

where R = e hw is the proper radial coordinate. Assuming the standard periodicity 2tt for the 
angular coordinate, there is generically a conical singularity at the center 

A — 27r ^1 - - J^) . (4.53) 

This conical singularity at R = can be seen as our brane world. The ordinary 3-space metric 
cannot depend on the internal space coordinate and therefore it is impossible to realize a warped 
structure. 

4.3 Cosmological behaviors on the brane 

Let us now discuss the cosmological behaviors of our solutions. In our real Universe, the scale 
factor is not oscillating. Thus, in this subsection we do not consider the solutions with oscillating 
scale factor and focus only on the sinh solutions of G w = G^ with f 2 < 0. Without loss of 
generality, we may assume f\ > 0. For simplicity, we discuss the behaviors in the large t limit: 



n « - I J 2(l - —gfj - 5-AgA fct =: p M 



^gfj - 1 + ^ ] ht =: p 1 f 1 t, 



4 



v = w « _I[ J 2(1 - ) + 1 + 4^ ) fit =: pjtf, (4.54) 
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Figure 1: The plot for p is shown as a Fi S ure 2: The P lots for Pi ( the solid curve ) 
function of g& for /i > and P v (^ ne dashed curve) are shown as a 

function of for /i > 0. 

where we also define := G^//i. Note that \g^\ < a/3/4 « 0.433. For -0.196 < g^ < 0.410, 
Pi > and for g^ > —0.395, p v < (see Fig. 1 and 2). 

In one choice of the proper time coordinate, dT = —e u °dt, the approximate spacetime metric 
is given by 

ds 2 =(- T) 2qv (dR 2 + R 2 (l- ^fde 2 ) -dT 2 +(- T) 2qi X>x*) 2 , (4.55) 

i=l 

where unimportant constants are eliminated by the proper rescalings of R and x l and A is 
defined in Eq. (4.53). We also defined the powers 

Qi ■= - 1 — r, Qv ■= - 1 — r- (4-56) 

bo I IpoI 

In the other choice of the proper time coordinate, dT = e u °dt, the approximate spacetime metric 
becomes 

A 3 
ds 2 = T 2Qv (dR 2 + R 2 (1- — ) 2 d0 2 ) - dT 2 + T 2gi J2( dxl ) 2 - ( 4 - 57 ) 

i=l 

In Fig 3, we show the behavior of q± and q v . For —0.196 < g<j, < 0.410, q\ < 0. For g^ > —0.395, 
q v > 0. For —a/3/4 < g^ < —0.335, qi > q v , while for the rest q\ < q v . 

For -0.196 < g<j> < 0.410, q 1 < and q v > 0. Thus, in the metric (4.55), in the T ->> 0- 
limit, the size of the ordinary 3-space increases and diverges within a finite time, while that of 
the internal space is shrinking to zero. 

For —a/3/4 < g^ < —0.335, q x > q v and q\ > 0. Thus, for this range of g^, the metric (4.57) 
describes the ordinary 3-space which is expanding faster than the internal space dimensions. 
In particular, for g^ < —0.395, the internal space is contracting. But since q\ < 0.406, an 
accelerating expansion, as in the inflationary or dark energy Universe, cannot happen. 
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-0.2 



Figure 3: The plots for q x (the solid curve) 
and q v (the dashed curve) are shown 
function of for f\ > 0. 

The Einstein frame: Let us briefly discuss how our solution behaves in the Einstein frame. 
Rewriting the original metric in this form, 

ds 2 = g$dx»dx v + e 2v dy 2 + e 2w d6 2 , (4.58) 

we obtain 

R = + 2v 2 + 2w 2 + ■ ■ ■ , (4.59) 

where R^ is the Ricci scalar associated with the metric g^. The gravity action can reduce to 

J d & x^GR ~ V 2 J d 4 x^g~^e v+w R {4) ~ V 2 J d 4 x^g^R (E) (4.60) 

where V 2 — J dvdw is the comoving volume of the internal space, which is assumed to be finite 
for instance by introducing a cutoff. The Einstein frame metric is given by 

3 

ds% := g [E) dx»dx v = e v+w g$dx»dx u = - e v+w+2u °dt 2 + e v+w+2ui ^(dx*) 2 . (4.61) 

i=i 

Using the solutions, we find 

1 4 

v + w + 2u! = -(v + w + 2u ) = -In sinh[/i(t - t x )] - -Gfb. (4.62) 
o o 

Thus, in the proper time coordinate system dr = ±e v+w+2u °dt, the Einstein frame metric can 
be written as 

3 

ds% « -dr 2 + ( =F r) 2/3 ^(cfe*) 2 , (4.63) 

i=l 

which corresponds to a contracting or expanding Universe filled by the stiff matter. 
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5 Conclusions 



We have investigated the dynamical 3-form flux compactifications and their implications for 
brane world cosmology in the six-dimensional Nishino-Salam-Sezgin model. We take the back- 
ground of the 3-form field acting on the internal space and timelike dimensions without the 
U(l) gauge field strength. 

The first class of solutions we discussed was the dynamical generalization [15] of the static 
solutions, recently obtained in [14]. It turned out that the dynamical generalization is possi- 
ble only for the special case. In this class, we found that the time evolution is restricted in 
the azimuthal dimension of the internal space. The ordinary three-dimensional space is not 
dynamical and hence does not give rise to a cosmological evolution. There is a curvature sin- 
gularity at the boundary of the internal space, which is more severe than the conical one. To 
construct the brane world model, following [14], we first removed the singular part from the 
original solution and then glued the remaining piece of the spacetime with an identical copy at 
the codimension-one boundary wrapping the axis of the rotational symmetry. This boundary 
may be regarded as our 3-brane world. Generically, because of the presence of the 3-form field 
the Lorentz symmetry in the ordinary four-dimensional spacetime is broken. In general, such a 
brane can be supported by the matter with the energy density p, the pressure in the ordinary 3- 
space p and that in the azimuthal dimension pg, which would be different from the pure tension. 
Since the inclusion of the time dependence reduces the number of the parameters,the brane can 
be embedded by the matter with pg = —p (but still p ^ — p). However, at the place where the 
Lorentz symmetry is recovered, the boundary brane world can be supported only by the pure 
tension, i.e., p = pg = —p. Note the solutions reported here and related implications for the 
brane world model were obtained under several assumptions. It is very interesting to explore 
the time dependent solutions and associated brane world models without these assumptions. 

The second class of models was obtained by exchanging the roles of the radial coordinate 
and the time coordinate from the first class one. At the center of the internal space, there 
is a conical singularity which may be interpreted as co dimension-two our 3-brane world. But 
it does not seem to be possible to realize a warped structure in the external dimensions, as 
expected in analogy with the previous case. Except for the oscillating dynamical solutions, the 
solutions led to the expanding or contracting ordinary 3-space, depending on the choice of time 
direction. Among the expanding solutions, there are decelerating and accelerating ones: In the 
latter solution, the scale factor in the ordinary 3-space diverges within a finite time. In the 
Einstein frame, however, the Universe always followed the expansion law of the one filled by the 
stiff matter, irrespective of the choice of parameters. 
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A Dynamical black 1-branes 

In this Appendix, we discuss another class of the dynamical solutions with the nonvanishing 
3-form field strength. The solution represents the black 1-brane extended along the y direction 
in the presence of a cosmological constant. We now take the following ansatz for the metric 

dsl = e 2u ^ (-dt 2 + dy 2 ) + e 2v ^ r \dr 2 + r 2 ^ 2 ), (A.l) 

the dilaton <fi = (f>{t, y, r) and the 3-form field strength 

H = E(t,y,r), r dt Adr Ady. (A.2) 

The (t, y) coordinates cover the worldvolume of the 1-brane, while the rest does the transverse 
spatial dimensions. The equations for the 3-form field are given by 

(e- 2 "+ 2u - 2 V£ r )- = (e- 2u+2v - 2 *r s E <r ) ty = (e- 2u+2v - 2 *r 3 E, r ), r = 0. (A.3) 

Defining V = 2u + 2v , the integration of Eq. (A.3) gives 



21 



where c is an integration constant. Then, the diagonal components of Einstein equations and 
the dilaton equation of motion are given by 

3 
— i 
r 



e -2u 



( — u ,yy ~ ^ u ,y v ,y — u + 2v + V — 2uv — uV + 2u 2 + Av 2 ) + e 2v ( — U ir V r — w rr w r ) 



-lull ^ C -V-2vjp , 9 2 



= -e~ 2u <f> 2 - He- v - 2v E, r + (A.5) 
2r 6 k 

e ~ 2U ( u ,yy ~ 2v ,yy ~ V ,yy + 2u ,y v ,y + u ,y V ,y ~ 2u2 y ~ Av l + u + u(V - 2v + 2ii)) 

+e~ 2v ( - u r V r - u rr - -u r ) = e" 2 > 2 - p-e~v-^ E + ^e*, ( A.6) 



-2u 



(v + v(V — 2u + 2v) — v m — v t y(V >y — 2u tV + 2v t y)j 



q 2 2 

+e~ 2v (-V rr -v rr --v r + v r V r - 2{u\ + v 2 T )) = e" 2 > 2 - ^ e - v ~ 2v E r + 9 —e\ (A.7) 
e~ 2u (v + v(V -2u+ 2v) - v m - v, y (V y - 2u, y + 2v jV j) 

+e-* ( _ V rr - l v _ ( w + i) Kr ) = ^-v-a^ + (A .8) 
v r r 2r 6 k z 



-2 a 



( - - (V - 2u + 2v)<p + + {V y - 2u y + 2v y )cf> >y ) 



+e~ 2v {<f>, rr + {V r + ~) 0, r ) - ^e"^ 2 ^ - = 0. (A.9) 

We require in Eq. (A.5) and (A.6) 



K 



2 



u r = -^—ce~ v E. (A.10) 
2r A 



Similarly in Eqs. (A. 8) and (A.9), 



k2 -v jp ± _ k2 „. 



v <r = -—ce- v E, cf) jr =—ce- v E, (A.ll) 

where we assume that V does not depend on r. Eliminating r-derivative terms in Eq. (A.7), 
by Eqs. (A. 10) and (A.ll), we obtain 

a 2 

e 2v-2u ^ + v(y-2u + 2v) - v m - v, y (V y - 2u y + 2v, y )] - ^ +2v 

L M~) (A.12) 



^3 y ' 

where we defined E := K 2 ce~ v E. Setting the right-hand side of Eq. (A.12) to be zero, we obtain 

E= ® r , H(t,y,r) = h(t,y) + %, (A.13) 
H(t,y,r) r z 

where Q is a constant. Then Eqs. (A. 10) and (A.ll) give that the metric functions and dilaton 
can be written as 

u = ~lnH{t,y,r), v = ^laH(t,y,r), = ~ In H(t,y, r), (A.14) 
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where the integration constants are set to be zero. It is straightforward to check that the 
off-diagonal components of the Einstein equation 



Au^ r v — u iT — 3v iT — ir = 0, 

4u, r V % y — U^ y — 3t\ ry — 4>,r4>,y = 0, (A. 15) 

are satisfied. We find that the time dependent parts of Eqs. (A. 5), (A. 6) and (A. 7) reduce to 

h yy + 3h = -^, (A.16) 

-h-3hyy=^r, (A. 17) 

4o 2 

h - hyy = -\, (A.18) 

respectively. Both of Eqs. (A. 8) and (A. 9) reduce to the same equation as Eq. (A.18). The 
solution of Eq. (A.16)-(A.18) is given by 

Q 2 

h{t, y) = ^ (t 2 - y 2 ) + a t t + a y y + a (A. 19) 

where a, (i = t, y, 0) are integration constants. Then, the left-hand side of the remaining 
off-diagonal component of the Einstein equation 

-iv t yV - 4v >y + Au tV v + Auv >y - <j><f) iy = 0, (A.20) 

is proportional to h iV , which vanishes because of Eq. (A. 19). As seen from Eq. (A. 19), h is 
no longer linear in the worldvolume coordinates but quadratic in them, which also happens in 
p-brane solutions with trivial or vanishing dilaton [22]. 
Let us summarize the solutions obtained here: 

u(t, y, r) = -v(t, y, r) = ~ In + ^(t 2 - y 2 ) + a t t + a y y + a ) , (A.21) 

cf>(t, y, r) = -\ In (% + 9 —{t 2 - y 2 ) + a t t + a y y + a ) . (A.22) 

Because of the dependence on the spatial worldvolume coordinate y as well as on the time 
coordinate t, the brane direction is not compact and this solution may not be suitable for 
constructing the brane world. One might imagine that there may be a solution with a quadratic 
order dependence on time in the case of the 2-form field (i.e. 0-brane). For the dilaton coupling 
parameters in the NSS model, such a dynamical solution does not exist. But it does for the 
other special coupling parameters. Our solution can be interpreted as one of the dynamical 
p-brane solutions, in a class of the theory with a cosmological constant. These solutions will 
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be discussed in detail in [23]. In particular, among these solutions, the 0-brane solution in the 
four- dimensional spacetime behaves as a black hole embedded into a FRW Universe filled by 
the matter with the equation of state w — — |, where w is the ratio of the pressure with the 
energy density. 
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